Relaxation of the Electron Spin in Quantum Dots Via One- and Two-Phonon 

Processes 
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We have studied direct and Raman processes of the decay of electron spin states in a quantum dot 
via radiation of phonons corresponding to elastic twists. Universal dependence of the spin relaxation 
rate on the strength and direction of the magnetic field has been obtained in terms of the electron 
gyromagnetic tensor and macroscopic elastic constants of the solid. 
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I. INTRODUCTION 

Electron spin relaxation in solids is related to such im- 
portant applied problems as electron spin resonance and 
quantum computing. Many interactions contribute to 
the relaxation time of the electron spin in a semiconduc- 
tor quantum dot. In principle, all of them can be elimi- 
nated (such as, e.g., impurities, nuclear spins, etc.) but 
the interaction with phonons cannot. Thus, spin-phonon 
interaction provides the most fundamental upper bound 
on the lifetime of electron spin states. The existing meth- 
ods of computing electron spin-phonon rates in semicon- 
ductors rely upon phenomenological models of spin-orbit 
interaction, see, e.g., Refs. [il3,lllillllllllli,|l3- 
These models contain unknown constants that must be 
obtained from experiment. Meantime, the spin-orbit cou- 
pling determines the difference of the electron g-tensor 
from the unit tensor. The question then arises whether 
the effect of the spin-orbit coupling on spin-phonon re- 
laxation can be expressed via the difference between the 
electron gyromagnetic tensor gap and the vacuum ten- 
sor gaSai3. Since gap can be measured independently, 
this would enable one to compare the computed relax- 
ation rates with experiment without any fitting parame- 
ters. In this paper we show that this can be done under 
certain reasonable simplifying assumptions. We obtain 
spin-phonon relaxation rates due to direct one-phonon 
processes and due to two-phonon Raman processes. 

Zeeman interaction of the electron with an external 
magnetic field, H, is given by the Hamiltonian 



-fJ-B dap SaHp , 



(1) 



where /x^ is the Bohr magneton and S = cr/2 is the 
dimensionless electron spin, with aa being Pauli ma- 
trices. One can choose the axes of the coordinate sys- 
tem along the principal axes of the tensor gafi- Then, 
5a/3 = ga^afs- Perturbation of Eq. iQJ by phonons has 
been studied in the past H, Q, |^ by writing all terms of 
the expansion of gaf3 on the strain tensor, Uaf3, permitted 
by symmetry. This gives spin-phonon interaction of the 
form Aap-ypUapSjBp with unknown coefficients Aapjp. 
To avoid this uncertainty we limit our consideration to 
local rotations generated by transverse phonons. The ar- 
gument for doing this is three-fold. Firstly, the rate of 



the transition accompanied by the emission or absorp- 
tion of a phonon is inversely proportional to the fifth 
power of the sound velocity [ij. Since the velocity of 
the transverse sound is always smaller than the velocity 
of the longitudinal sound 1^, the transverse phonons 
must dominate the transitions. Secondly, we notice that 
interaction of the electron spin with a local elastic twist 
generated by a transverse phonon does not contain any 
unknown constants. Consequently, it gives parameter- 
free lower bound on the spin relaxation. Thirdly, for a 
dot that is sufficiently rigid to permit only tiny local rota- 
tions as a whole under an arbitrary elastic deformation, 
the emission or absorbtion of a quantum of the elastic 
twist is the only spin-phonon relaxation channel. 

The angle of the local rotation of the crystal lattice in 
the presence of the deformation, u(r), is given by 



-V X u 

2 



(2) 



As stated above, the gyromagnetic tensor gap is deter- 
mined by the local environment of the quantum dot. In 
the presence of long wave deformations of the lattice, the 
whole environment is rotated so that the gyromagnetic 
tensor becomes 

gap = ^aa'^pP'ga'P' , (3) 

where Wap is the 3x3 rotation matrix given by S4){r). 
One can thus write the total Hamiltonian in the form 



'H = -H'y + n 



ph : 



where 



- -l^B9a'P' i^aL^a) {^p'^pHp 



(4) 

(5) 



and TLph describes free harmonic phonons. In the above 
formulae, 5(p must be understood as an operator. Indeed, 
from the canonical quantization of phonons and Eq. Q 

one obtains 



(50 = 



8pV 



E 

kA 



[ik X GkA] e"' 



akA 



^ kA 



(6) 



where p is the mass density, V is the volume of the crys- 
tal, GkA are unit polarization vectors, A = ti,t2,l denotes 
polarization, and ujk\ = v\k is the phonon frequency. 
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II. DIRECT PROCESSES 

In order to account for spin transitions accompanied 
by the emission or absorption of one phonon one needs 
to consider terms up to first order in phonon amphtudes. 
Therefore, with the help of the expansion of the rotation 
matrix to the first order in Sep, = S^p — eap~f5(j)^^ we 
obtain the full Hamiltonian 

H ^Ho + Ti-s-ph , 'Ho ~ Jiz + 'Hpu , (7) 



where 



(8) 



Spin-phonon transitions occur between the eigenstates of 
Tio, which are direct products of spin and phonon states 



(9) 



Here |V'±) are the eigenstates of Hz with energies E± 
and \ 4>±) are the eigenstates of Hph with energies Eph±. 
Spin-phonon processes conserve energy, _E+ + -Eph+ ~ 
+ Ep\,^. For Hs-ph of Eq. (Q, the states \4>±) differ 
by one emitted or absorbed phonon with wave vector k. 
We will use the following designations 



(10) 



To obtain the relaxation rate one can use the Fermi 
golden rule. The matrix element corresponding to the 
decay of the spin |^'+) — > l^*-) can be evaluated with 
the help of Eqs. ijHl and JB)) we obtain 



{^Ans-ph\^+) 



kA 



ViA 



: K • [k X GkA] 



where the components of vector K are given by 



ill) 



(12) 



Note that only the transverse phonons are considered in 
the summations. For the direct process, the decay rate, 
Td, is, then. 



Td= Focoth 



HujQ 
2fcsT 



where 



huj(] 



E^ 



E- 




is the distance between the two spin levels, and 
To = ^ |VkAp27rJ(tJkA - (^o) ■ 



(13) 



(14) 



(15) 



Using Eq. ifTT)) and replacing ^j^. by {V/{2ny) J (Pk, one 
obtains 



To 



1 IKp^g 
12nh pvf 



(16) 



where Vt is the velocity of the transverse sound. A 
straightforward calculation of the spin matrix elements 
yields 



IKI- 



^ ~ 



(17) 



Then, the decay rate can be written in the final form 



3ttp \ hvt J 



Erin) , 



where n = H/H and 



3/2 



coth 



E 



2 2 



a0 



nl + n} 



(18) 



1/2- 
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2^2 ^2 



(19) 



kA 



If the field is directed along the Z-axis, Eq. (|19|) simplifies 
to 

FAe.) = |[(ff. - g.f + {g. - gy?] coth {^^) • 

(20) 

When components of gap are of order unity and ksT > 
PgH, then Tu ^ {kBT/h){figH/Et)'^, where we have 
introduced Et = {h^pvfY^'^ - lO^K. 



III. RAMAN PROCESSES 

Spin-lattice relaxation by Raman scattering is a two 
phonon mechanism consisting of a spin transition accom- 
panied by the absorption of a phonon and the emission of 
another phonon of different frequency. In spite of being 
a second order process, its contribution can be very im- 
portant, since the phase space of the phonons triggering 
the transition is not limited to the distance between the 
levels as in the direct case. To describe such processes, 
we need to consider terms up to second order in phonon 
amplitudes in the Hamiltonian. To this end we expand 
the 3x3 rotation matrix to second order in (50 

= Saf3 - eal3y6(j)^ + - [5(t>^6(j)p - (5q/5 ((50)^] (21) 

and obtain the full Hamiltonian from Eqs. Q and Q 
'H = 'Ho+ T-is-ph , T-is-ph = H^^lph + n'^s-ph ' (22) 
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where Ti-l^lpf^ is given by Eq. 



and 



n 



(2) 



(23) 



Raman processes may dominate over direct processes 
for a small energy difference between the spin states, 
hiuo ^ ksT. Thus we will consider terms to the lowest 
order in H. For simplicity we will study the case where 
the field is directed along the Z-axis. Then the matrix 
element becomes, 



Again, we will study the spin-phonon transitions be- 
tween the eigenstates of Hq, |^'±) — |^±) §5 |0±)- Here, 
again, I ^/;^) are the eigenstates of Hz with energies E± 

and the eigenstates oiHph with energies Eph±- We 
consider Raman processes, in which a phonon with wave 
vector k is absorbed and a phonon with a wave vector q 
is emitted. We will use the following designations for the 
phonon states 



Mr, 



(2) _ Mjjg 



M 



[{2gy - 5, - g,)M^^_j, 

ph-Ri ' (27) 



where Af"/' ^ = M""^ 



'^ph-R 



ph-R 



M: 



pa 

ph-R 



with 



_) = |nk,71q), = |71k-l,?lq+l). (24) = (riq + 1 1 5 0„ | ^q) ( rifc - 1 | (5<?!>0 | rifc ) 



To obtain the relaxation rate of the transition |^'+) 
\^ -) one needs to evaluate the matrix element of the 
process, which is the sum of the matrix element with 

in the second order: M^i = 



^2 [k X ekAj„ [q X GqAj 



8pV 



fif^,. and that with Ti'^ 



(riq + l)nk 
(28) 



M 



and 



>—ph 
(2) , 



M 



(1+1) 



where 



M 



(1+1) 



E 



M 



(2) 



s—ph 



(vl/_|7ifj^j*+) 



(25) 



The Raman rate of the transition l^*-) can be 

obtained by using the Fermi golden rule. 



K 



(1) 

S — ph 



(1) 



-ph 



Ti? = [{2gy - - gz)"^ + {2gx - gy - gz)"^] 



3024 h \ Et J 



Et 



(29) 



(1) 

S—ph 



E+ 



(1) 

S—ph 



(26) 



The intermediate phonon states are |nk — 1,'^q) in the 
first term and |nk, Wq + 1) in the second term. 



The ratio of Raman and direct rates is of order 
Th/Td ~ 10-^ {kBT/iiBH)^{kBT/Ety. Consequently, 
Raman processes dominate over direct processes at high 
temperature and low field. 

This work has been supported by the NSF Grant No. 
EIA-0310517. 
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